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AD3TRACT 

For  a  cylindrically  symmetric  plasma  column  whose  electron 
density  N  is  a  slowly  varying,  monotoni tally  decreasing  function  of  radius 
r,  it  is  possible  to  calculate  explicitly  both  the  diffraction  pattern  from  the 
knowledge  of  Nfr)  and,  conversely,  N(r)  from  the  knoudedge  of  the  diffrac¬ 
tion  pattern.  If  the  diffraction  pattern  is  obtained  experimentally,  N(r)  can 
be  calculated  by  a  cumbersome  numerical  procedure.  Instead  of  doing  this, 
the  diffraction  pattern  can  be  approximated  by  one  of  a  family  of  convenient 
analytical  expressions  for  which  the  integration  can  be  carried  out  easily. 
Alternatively,  one  can  attempt  to  infer  N(r)  by  assuming  a  functional  form 
for  N(r)  with  one  or  more  parameters,  calculate  the  diffraction  pattern  and 
compare  it  with  the  observed  one. 
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1.  Introduction 

It  i#  |)Ossible  to  determine  the  electron  densi'y  distribution  in  a 
lossless  plasma  cylinder  by  illuminating  it  ^*ith  a  plane  electromagnetic 
wave  and  observing  the  diffraction  pattern.  Let  us  assume  that  the  plasma 
cylinder  is  long  compared  to  wavelength  and  that  the  plasma  density  de[  ?nds 
on  the  distance  from  the  axis  alone.  This  plasma  cylinder  is  illuminated  by 
a  weak  plane  electromagnetic  wave  incident  at  right  angler  to  the  axis, 
polarized  so  that  its  electric  vector  is  parallel  to  the  axis.  The  electric 
field  has  a  single  axial  component  which  satisfies  the  wave  equation 


whe  re 


and 


[-7^  ♦  k^r)l  =  0 


2  2  , 

=  k  -N|r)u  e/m 
o  o 


(1) 

-f' 

(2) 


w  -  angular  frequency  of  the  wave 

^  •  c  •  magnetic  permeability  and  dielectric 
o  o 

constant  of  free  space 

e,  m  •  electronic  charge,  mass 

N(r)  •  electron  density  (no.  of  electrons/m^) 

k  -  w/  p  (  .  free  space  wavenumber 

o  y  o  o 


Far  away  from  the  plasma  cylinder  the  field  consists  of  the  incident  plane 
wave  and  an  outgoing  cylindrical  wave.  From  the  dependence  of  th^  power 
density  in  the  cylindrical  wave  on  the  angular  direction  or  on  the  frequency 
It  is  possible  to  deduce  the  radial  distribution  N(r).  This  inverse  scattering 
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problem  has  received  considerable  attention  in  quantum  mechanics.  The 
wave  equation  then  is 


whe  re 


2 

k  ir)] 


4,-0 


2 

k  (rj  = 


2m 


2  m 


V(r) 


The  energy  E  here  plays  the  same  role  aa  the  frequency  does  in  the  case 
of  plasma,  and  the  potential  energy  V(r)  corresponds  to  the  electron  density 
distribution.  In  quantum  mechanics  one  normally  tries  to  construct  the 
potential  function  V(r)  from  known  scattering  data  (either  total  scattering 
cross-section,  or  differential  scattering  cross -section  in  a  fixed  direction) 
as  a  function  of  energy.  This  is  basically  the  same  procedure  as  is  used, 
for  example,  in  ionospheric  sounding,  where  frequency  is  varied.  The 
method  dealt  with  here  utilizes  scattering  data  in  all  directions  at  a  single 
energy  (frequency). 

Whatever  type  of  data  is  used,  the  exact  inversion  problem  is  an 
extremely  diffirull  one  to  solve^.  For  this  reason  we  may  resort  to  an  ap¬ 
proximation.  U  the  ‘'ariation  of  properties  of  the  medium  is  slow,  i.e.,  the 
relative  change  per  wavelength  is  small,  then  we  can  resort  to  a  WKB  or 
geometric -optical  approximation.  We  shall  assume  that  the  geometric - 
optical  approximation  is  valid,  calculate  the  differential  scattering  cross- 
section  from  the  electron  density  distribution  on  that  basis,  and  then  invert 
that  relationship  so  as  to  obtain  the  electron  density  distribution  from  the 
differential  scattering  cross-section. 


1.  I.M.  Celfand  and  B.  M.  Levitan,  lav.  Ahad.Nank  SSSR,  Mat.  Ser.  IS,  309 

(1951) 

Ooktady  Ahad.  Nank  SSSR,  77,  S57(1951) 
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2.  Differencial  Scattering  Croae  -  jection  and  Impact  Paran'.eter 

Consider  a  plane  wave  incident  perpendicularly  on  a  plasma  column 
the  rays  are  initially  parallel  straight  lines  (See  Fig.  la).  When  these  rays 
enter  the  ionized  region,  they  are  deflected.  The  im';ict  parameter  b  for  a 
given  ray  is  defined  as  the  distance  between  the  central  ray  (ray  in  the  plane 
of  symmetry  of  the  plasma  column)  and  the  given  ray.  If  the  incident  wave 
carries  unit  power  per  unit  width  (this  is  a  '  -o-dimensional  geometry),  then 
the  power  contained  between  rays  with  impact  parameter  a  and  b  t  db  is 
db.  If  these  rays  emerge  from  the  plasma  at  angles  Q  and  6  +  d0,  then, 
on  the  basis  of  geomet  ric -optical  consideration, 

<rd0  i  db  (4) 


where  9  is  the  scattered  power  per  unit  a.ngle,  or  the  differential  scattering 
cross-section.  We  see  that  j(6)  and  0(b)  are  simply  related  as  follows: 


db 


l‘>) 


The  firnt  step  in  the  inversion  process  is  to  calculate  0(b),  given  (r(0). 

Even  if,  as  we  shall  assurr>e,  the  ionization  density  N  is  a  monotonic  de¬ 
creasing  (unction  of  r,  there  may  be  a  difficulty  in  this  step.  We  have  to 
distinguish  between  soft,  or  penetrable,  obstacles,  and  hard,  or  impenetrable 
ones.  If  the  plasma  column  is  impenetrable  then  the  central  ray  is  deflected 
through  180°,  while  distant  rays  are  not  deflected  (see  Fig.  lb).  In  that 
case,  b(0)  is  a  single -valued  function,  and  there  is  no  difficulty  in  inte¬ 
grating  e(9)  to  obtain  b(0),  and  then  to  invert  that  relation  to  obtain  6(b). 

If,  on  the  other  hand,  the  plasma  column  ie  penetrable  (Fig.  Ic),  the  central 
ray  is  not  deflected,  and  other  rays  are  deflected  through  various  angles 
ranging  from  ssro  to  a  certain  riiaximum  value  9  (maximum  diffraction 
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angle).  Hence  there  are  two  ray»  deflectec  in  any  given  direction  9,  and 
b(9)  ' 3  double -valued. 

In  the  latter  case  the  differential  scattering  cross -sectioit  will 
oscillate  rapidly  as  a  function  of  9  as  a  result  of  interference  oetween 
the  two  rays  deflected  in  the  same  direction.  At  the  peaks  the  two  intensities 
add,  at  the  minima  they  subtract.  If  we  plot  intensity  (square  root  of  dif¬ 
ferential  cross-section)  vs.  9.  we  can  readily  construct  the  curves  showing 
the  variation  of  the  sum,  I  ,  a.nJ  the  difference,  I  ,  of  the  two  ray  inten¬ 
sities  Ij  and  (see  Fig.  2).  Assuming  >1,,  we  have 


-  I  ) 


1  t 

j(I  +  I  ) 


(6) 


There  are  difficulties  to  be  expected  in  reconciling  experimental  data  with 

ray  theory.  If  we  start  from  a  dependence  of  the  impact  parameter  b  on  9^ 

such  as  shown  on  Fig.  2c  then  we  notice  that  the  intensity  associated  with  the 

inner  ray  (1)  becomes  infinite  at  9  ,  and  the  intensity  associated  with  the 

o 

outer  ray  (2)  becomes  infinite  at  9  =  0  and  at  0  .  No  such  singularities 
will,  of  course,  be  observed.  They  appear  in  the  geometric-optical  approxi¬ 
mation  because  the  approximation  is  not  valid  at  tho^e  values  of  9  as  a 
consequence  of  the  behavior  of  the  caustic  surface  in  this  case.  Should  the 
ionization  density  be  identically  zero  outside  some  finite  radius  then  the 
singularity  at  9  ^  0  may  disappear.  In  any  case,  we  can  utilize  only  those 

prrtions  of  the  curves  which  are  not  too  close  to  9  =  0  or  9=9.  The 

o 

problem  is  then  to  construct  all  of  the  curve  on  Fig.  2c  from  the  knowledge 
of  the  slope  of  this  curve  (or,  rather,  the  slopes  of  its  two  branches)  in 
the  range  a  <  9  <  9^  - 
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It  can  be  shown  that,  for  small  values  jf  impact  parameter,  the 

ray  deflection  is  proportional  to  the  impact  parameter.  This  means  that  we 

can  extend  the  curve  to  9  =  0  simply  by  assuming  that  remains 

constant  for  0  <  9  <  a.  Thus  we  obtain  b  vs,  9  from  the  origin  to  9  =  9  - 

o 

We  shall  assume  that  while  the  scattering  data  in  the  range  9  -  P  <  9  <  9 

o  o 

cannot  be  used,  the  value  of  9^  is  known,  and  that  the  value  impact  param¬ 
eter  at  9  =  9q  -  P ,  as  well  as  the  slope  of  the  b  vs.  9  curve  (the  slope  is 
simply  the  differential  scattering  cross-section  a),  are  known,  and  finally 
that  the  shape  of  the  9  vs.  b  curve  is  parabolic  (see  Fig.  2d),  i.e., 

(9  -  9)  «:  (b  -  b)^ 

o  o 

With  these  assumptions,  the  difference  between  the  two  values  of  impact 

parameter  at  9  =  9  -  p,  b  -  b  is  given  by 
o  2  1 

b^-  bj  =  4p(r(9^-  P)  (7) 

where  (r(9  -  P)  is  the  differential  scattering  cross-section  in  the  direction 
o 

9  =  9  -  p.  Since  b  is  known,  equation  (7)  establishes  a  point  on  the 
o  i 

second  branch  of  the  b  vs.  9  curve.  From  that  point  we  can  continue  the 
curve  by  using  experimental  values  of  scattering  cross-section  again. 

3.  Ray  Path  and  Deflection  Angle 

The  differential  equation  satisfied  by  the  ray  path  r(9)  in  a 
cylindrically  stratified  medium  can  be  readily  obtained  from  Snell's  law. 

If  the  dielectric  constant  is  a  function  of  r,  then  Snell's  law  takes  on  the 
form 

r  sin  a  =  constant  (8) 

where  a  is  the  an^le  between  the  ray  direction  and  the  radial  direction 
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{see  Fig.  3).  Since  for  large  r,  t— rsina— ♦b,  the  constant  on  the 
right-hand  side  of  (8)  is  simply  the  impact  parameter  b  of  the  given  ray. 
From  triangle  PP'Q  (see  Fig.  3)  we  see  that 

2 

=  tan  a  .  (9) 

Eliminating  a  between  equations  (8)  and  (9),  we  obtain  the  differential  equa¬ 
tion  of  the  ray 

(rd9/dr)^  =  (r^f/b^c  -  1)'^  (10) 

o 

In  the  case  of  a  lossless  plasma  the  dielectric  constant  is 


c/e 

o 


2  2,2 

ik/k  )  =  I  -  N(r)e  /c  mu 

o  o 


(H) 


Let  us  introduce  a  notation  which  has  the  double  purpose  of  making  the  rela¬ 
tion  (ll)^and  those  which  folloW|more  compact,  and  at  the  same  time  bring 
out  the  similarity  of  the  problem  at  hand  and  the  corresponding  one  in  particle 
mechanics;  let  us  introduce 


V(r) 


N(  r)e 


2  .  2 
f  E  TTlw 
O 


(12) 


Combining  equations  (10),  til)  and  (12)  yields  the  final  form  of  the  differen¬ 
tial  equation  of  the  ray  with  the  impact  parameter  b 


|S£|  . 

The  angle  between  the  direction  along  which  the  ray  comes  in  from  infinity, 
and  the  direction  along  which  it  goes  out  to  infinity  again  is  tlien 


dr 


dr  . 
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The  deflection  angle  of  the  ray,  i.e.,  the  limiting  value  of  the  angular  co¬ 
ordinate  as  r— »aa,  becomes 

?  -2  r  -2  -2  -2  1 

0(b)  =  w  -  Z  J  r  b  -r  -b  V{r)  dr  (14) 

*^o 

where  r  ,  the  distance  of  closest  approach  for  this  ray,  is  the  zero  of  the 
°  dr 

expression  in  brackets  (-—  =  0  at  r  ). 

‘d0  o 

We  will  now  investigate  the  function  6(b)  ^or  various  ionization 
density  distributions.  We  have  to  distinguish  between  two  types  of  distribu¬ 
tions;  penetrable  and  impenetrable.  The  illustrative  examples  will  be  of 
the  former  type.  We  shall  consider  three  ionization  distributions; 

(1)  uniform  (2)  parabolic  and  (3)  modified  paraboUc.  The  mathematical 
definitions  of  these  three  distributions,  as  well  as  the  relations  between  the 
deflection  angle  and  the  corresponding  impact  parameter  are  .'isted  in 
Table  1.  The  functions  9(b)  were  obtained  simply  by  evaluating  the  integral 
in  Eq.  (14)  with  the  appropriate  distribution  function  V(r).  The  results 
listed  in  Table  I  are  also  shown  graphically  on  Figs.  4,  5,  and  6.  The  case 
of  the  uniform  cylinder  is  different  in  character  from  the  two  others  because 
the  medium  is  discontinuous  in  this  case.  The  geometric  optical  field  con¬ 
sists  of  refracted,  reflected  and  many  multiply  reflected  and  refracted  con¬ 
tributions.  If  we  consider  only  the  refracted  ray,  then  the  deflection  angle 
increases  monotonically  until  b  s  a/1  -  ;  for  larger  values  of  impact 

parameter  the  incident  ray  is  totally  reflected,  until,  for  b  >  a,  the  in¬ 
cident  ray  completely  misses  the  cylinder  and  is  not  deflected.  The  remain¬ 
ing  two  models  present  no  such  complications.  In  the  case  of  a  parabolic 
distribution  the  deflection  angle  is  identically  zero  if  the  ray  misses  the 
cylinder,  i.e.,  b  >  a,  while  in  the  modified  parabolic  distribuMon  in  which 
the  extent  of  plasma  is  infinite,  the  deflection  angle  tends  to  zero  as  b  in¬ 


creases. 
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TABLE  I 
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It  is  interesting  to  observe  the  dependence  of  the  maximum  ray 

deflection  angle  on  maximum  density  V  for  the  two  continuous  models. 

o 

The  maximum  deflection  angle  should  be  easily  observable  since  outside 
this  angle  there  is  no  geometric  optical  scattering,  only  true  diffraction. 

At  sufficiently  high  frequencies,  the  differential  scattering  cross-section 
should  drop  rapidly  to  very  small  values  at  this  angle.  We  see  from  Fig.  7 
that  there  is  very  little  difference  between  the  two  models  in  their  effect  on 
the  maximum  deflection  angle.  Therefore  we  may  conjecture  that,  so  long 
as  the  variation  of  density  near  the  center  of  the  plasma  column  in  quadratic, 
the  maximum  density  can  be  inferred  from  the  maximum  deflection  angle 
with  reasonable  accuracy.  Since  the  relation  between  the  two  quantities 
is  particularly  simple  in  the  case  of  the  parabolic  model,  we  can  make  use 
of  it,  thus  concluding  that 

V  sine  (15) 

o  o 

may  be  a  fair  approximation  for  a  great  variety  of  smoothly  varying  distribu¬ 
tions. 

4.  The  Inverse  Problem 

We  shall  now  proceed  to  invert  equation  (14),  i.e.,  to  solve  it  for 
V(r)  if  e(b)  is  given.  In  order  to  accomplish  this,  we  shall  reduce  the 
equation  to  one  of  Abel  type,  and  then  use  the  well-known  solution  of  the  Abel 
integral  equation.  The  procedure  used  here  is  identical  with  that  employed 
in  a  previous  solution  of  the  integral  equation  (14).^  It  is  convenient  to  intro¬ 
duce  the  new  variables  x,  u,  v,  and  w,  defined  as  follows: 


2.  J.B.  Keller,  I.  Kay  and  J.  Shmoys,  'Determination  of  the  Potential  from 
Scattering  Data*  ,  Phys.  Rev.,  Vol.  102,  No.  2,  557-559,  April  15,  1956. 
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X  = 

b-^ 

(16a) 

-1 

u  = 

r 

(16b) 

v(u)  = 

1  -V 

(16c) 

2  -1 

W  3 

U  V 

(I6d) 

In  ternTiS  of  the  new  variables,  equation  (14)  becomes 


ir  -9(x) 

2 


/ 


V  ^^^(du/dw) 


(x-w) 


1/2 


dw 


(17) 


If  we  now  recognize  that  the  numerator  of  the  integrand  is  simply  an  unknown 
function  of  w,  equation  (17)  is  of  Abel  type.  Denoting  the  numerator  of  the 
integrand  in  (17)  by  g(w),  we  have 


V  *^^(du/dw)  =  g{w)  =  ^ 


w 


y  -  e(x) 

2w  •'  ,  TJl 

o  (w-x) 


hs- 


dx 


w 


1  -1/2  1  d 

S  —  W  -  - 

2 


e(x)dx 


2»  dw  J  l/z 

(w-x) 


From  (16d)  we  see  that 


du 

dw 


1  -1/2  1/2  I  1/2  -1/2  dv 

-w  '  V  +  -w  V  — 

2  2  dw 


and  therefore 


g(w)  »  J 


—  w 


(18) 


(19) 


Solving  this  equation  for  v  in  terms  of  g(w)  and  noting  that  v  a  1  when 
w  >  0,  we  obtain 
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w 


'i  o/w'I  -  <w'I  1  dw' 


v(w)  =  exp  j  2(w')  g(w')  -  (w')  dw'  .  (20) 


Substituting  the  expression  for  g(w)  from  (18)  into  (20),  we  obtain 


v(w)  =  exp 


w 

-1  r  1  g  f  3(x;ax 

»  •'  r—  dw*  ■'  f'  ;■  — 
o  /w'  o  /w'-x 


(21) 


d  t  8(x)dx 
o 

Having  thus  found  v(w),  we  can  now  calculate  r  =  u"^  ”  v'^/^w“  1/2  and 
V  =  1  -  V.  The  calculation  of  V(r)  is  then  a  parametric  one,  with  w  as  a 
parameter. 


5-  Spe  cial  Cases  of  9  vs.  b  Variation 

In  principle  the  procedure  of  the  preceding  section  can  be  ap> 
plied  to  numerical  data  for  0(b).  The  integral  in  (21)  would  have  to  be 
evaluated  numerically  for  a  large  number  of  values  of  parameter  w; 
this  is  clearly  a  cumbersome  procedure.  If,  instead  of  evaluating  the  inte¬ 
gral  numerically,  we  could  approximate  the  numerical  data  for  9(b)  by 
one  of  a  family  of  analytic  expressions  which  would  permit  explicit  evalua¬ 
tion  of  v(w),  the  procedure  would  be  greatly  simplified.  We  shall  now  dis¬ 
cuss  one  such  family  of  functions  9(b).  These  functions  are  linear  for 

small  value  of  b,  reach  a  maximum  9  =  9  at  b  =  b  (this  implies  a  soft, 

o  o 

or  penetrable  distribution),  and  then  tend  to  zero  as  some  inverse  power  of 
b.  The  defining  equation  for  9  is 


9=0  — 
o  b 

o 


i  ♦  n 


[n  ♦  {W\)  J 


nfl 

2 


(22) 


For  large  b,  9-^b  We  are  going  to  consider  the  special  cases  n  «  1,2. 
(See  Fig.  8.) 
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For  the  case  n  =  1  we  have 


So  that 


J_  _  2b/bo 

®o  I  +  (b/b  ^ 
o 


e(x)dx 


o  y w-x 


1  +  w  b 


v(w)  =  exp  -0  b  f  (w*  +  b  )  (W)  dw' 

o  o  ^  o 

o 

=  exp  1^-20^  /(wb^  +  1 )]  ^  j 


The  parametric  equations  for  the  normalized  plasma  density  V  in  this  case 


are  then 


V(w)  =  1  - c 


,  2-1/2 
-20  (1  +  1/wb  ) 
o  o 


.  .  -1/2  0o(l  +  l/wbo) 

r(w)  »  w  c  ” 


It  is  clear  from  the  above  equation  that  small  values  of  w  correspond  to 
large  values  of  r,  and  conversely.  While  it  i.t  impossible  to  eliminate  w 
between  equations  (26)  and  (27),  and  thus  obtain  V  as  a  function  of  r 
directly,  it  is  possible  to  accomplish  this  approximately  for  very  small  or 
very  large  r.  We  are  primarily  interested  in  the  value  of  V  at  r  a  0, 
i.e.,  V^,  and  in  the  manner  in  which  V  tends  to  zero  for  large  r; 

w-«ao,  r— #0,  (28) 
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1  'Z  /wb  ^ 

w— ►O.r^jw  ,Vsil-£  °  °;s29  J\uh  s;  29  b  /r 

u  o  o  o 


For  the  case  n  =  2  we  have 


(29) 


9  b 
o  o 


2  +  (b/b^)  J 


3/2 


(30) 


so  that 


ard 


,3/2^ 

w  3  9 

9(x)dx  o 


/ 


O  V 


2b 


1  -  2b  w  2y2b  w 
»  .  -1  o  o 

^  -  sin  ^  ^ 

I  +  2b  w  1  +  2b  w 
o  o 


v(w)  =  exp 


s  exp 


b'  w  ,  , 

— — /  (If  2w*  b  )  dw' 
If  ^  o 

Q 


ffC  +  l/2b^w) 
o 


(31) 


(32) 


The  parametric  equations  for  V(r)  in  this  case  are; 


V(w)  =  I  -  exp 


1/2 

r  w)»w  exp 


/ir)  (1  ♦  l/Zb  w) 
o 


■‘] 


2w)  (I  +  l/2b^w) 
o 


(33) 

(34) 


The  behavior  of  V(r)  at  r  s  0  and  r  •>«cs  is,  in  this  case 


w  — ♦  oo:  V 


.  i . 


0 :  r 


(35> 

z.  b' 

w3 - j-2-2-  (it) 

*  r 
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The  plasma  density  distribution  obtained  from  these  models  for  two  values 
of  maximum  deflection  angle  is  plotted  in  Fig.  9. 

We  can  use  the  results  obtained  in  this  section  to  study  further 
the  relation  between  maximum  deflection  angle  and  maximum  normalized 
plasma  density.  The  results  obtained  for  the  two  models  used  in  this  section 
are  compared  with  the  results  for  a  parabolic  distribution  in  Fig.  10.  The 
relation  between  the  maximum  deflection  angle  and  maximum  electron  density 
for  either  of  the  two  models  investigated  in  this  section  is  not  very  different 
from  that  obtained  for  a  quadratic  density  variation.  Given  the  frequency 
of  the  probing  wave  and  the  maximum  deflection  angle,  we  can  again  estimate 
the  maximum  electron  density  in  the  ionized  column  by  the  use  of  eq.  (15). 

If  more  accuracy,  or  more  information  about  the  distribution  is  required, 
we  must  compare  the  experime>ttally  observed  relation  between  a  scattering 
angle  and  impact  parameter  with  those  used  in  this  report,  pick  out  one  that 
corresponds  mere  closely  than  the  others,  and  look  up  the  corresponding 
electron  distribution. 
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